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Abstract 

Heavy flavor research is a vigorous and active topic in high-energy QCD physics. Comparing 
theoretical predictions to data as a function of flavor provides a unique opportunity to tease 
out properties of quark-gluon plasma. We explicitly demonstrate this utility with energy loss 
predictions based on the assumption of 1) a weakly-coupled plasma weakly coupled to a high- 
er probe using pQCD and 2) a strongly-coupled plasma strongly coupled to a high-p T probe 
using AdS/CFT; we find that while the former enjoys broad qualitative agreement with data, it is 
difficult to reconcile the latter with experimental measurements. 



1. Introduction 

Our goal as nuclear physicists is to quantitatively extract experimentally and understand the- 
oretically the properties of nuclear matter; as part of this goal, we wish to compose the phase 
diagram of the strong force. This is an extremely immodest goal. For example, detailed first- 
principles calculations for the phase diagram of hydrogen, the most simple QED system, is con- 
temporary research [1]. Nevertheless quark-gluon plasma (QGP) provides us with a unique op- 
portunity to probe experimentally and theoretically the emergent, many-body physics of a non- 
Abelian gauge theory in a certain region of its phase diagram. We have a number of tools at 
our disposal for exploring the properties of QGP from experimental measurements, for example: 
low-p T particles, electromagnetic probes, quarkonia, high-p T light hadrons, and high-p T heavy 
hadrons. Heavy flavor measurements and theory are an important piece in the puzzle we are try- 
ing to put together to form a consistent and coherent picture of heavy ion collision phenomena. 

Unfortunately length limits the breadth of the coverage here of the fascinating and impor- 
tant aspects of heavy flavor research in high-energy QCD; we will focus on the physics we can 
learn from high momentum particles. High-p T particles are especially interesting as they are the 
decay products of high-p T partons, which are the most direct probe of the relevant degrees of 
freedom in a quark-gluon plasma (QGP) [2, 3]. One is able in principle to learn about QGP by 
making an assumption regarding the physics of the QGP and comparing the necessary theoreti- 
cal consequences of those assumptions to data. One hopes to use this approach to falsify certain 
assumed descriptions of the plasma and add evidence for others. Just as one requires consistency 
of the entire picture of all data, one requires consistency of description of measurements associ- 
ated with energy loss. It turns out that this consistency is quite hard to achieve, and, as a result, 
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energy loss provides us with a valuable window through which to actually measure the physics 
of quark-gluon plasma. To explicitly demonstrate the power of this method we will take two 
extreme, generic assumptions regarding the medium and how it couples to high-p T probes and 
compare the results to data. 

Heavy flavor is especially interesting because of the additional experimental constraints it 
puts on energy loss calculations, and thus on the potential properties of quark-gluon plasma. In 
particular, an energy loss calculation, in a broad sense, results in a probability of a parton losing 
some of its initial momentum, P(Ap T \p T , L, T, Mq, R), where L is the pathlength the parent 
parton travels, T is the temperature of the plasma, Mq is the mass (or effective mass) or the parent 
parton, and R is the representation (i.e. is the parent parton a gluon or quark). Unfortunately one 
cannot alter these parameters experimentally to test energy loss theories; rather one changes, for 
instance, the collision species, the ^s of the collision, or — as is especially useful for investigating 
heavy flavor — the mass of the measured hadron, and compares to theoretical predictions. 

Qualitatively, one expects a simple ordering of the energy losses as one changes parent parton 
species, specifically 1 

AE b <AE c <AE u ^ s <AE g . (1) 

However, it is important to emphasize that energy loss ordering does not necessitate an or- 
dering of the nuclear suppression factors, Raa(Pt), because Raa involves a convolution over 
the fragmentation functions and, most important, the production spectrum. For approximately 
power law production of high-p T partons, as predicted by pQCD, and fractional momentum loss 
c = (Ptj - PT,f)/PT,i, one can show that [5] 

where the angular brackets refer to averaging over the geometry. As one can see from Eq. (2), 
the softer the spectra the larger is n and thus the smaller Raa for the same Apj . As seen in Fig. 1, 
the relative changes in production power laws nq for different parton species lead to a crossing 
in Raa for different hadron species (as first noted in [6]). 

The use of heavy flavor allows for an additional test of energy loss theory, and it turns out 
that by comparing to the centrality, center of mass energy, and especially flavor dependence of 
data yields a stringent constraint on the theory; see how, for instance, a simultaneous description 
of pion and non-photonic electron suppression at RHIC precludes the possibility of pQCD-based 
radiative only energy loss [7]. 

2. Extracting Physics by Comparing to Data 

2.1. Strongly-coupled Medium Strongly Coupled to a Probe 

To begin our energy loss comparison to data using two extreme assumptions about the physics 
of the QGP, let's consider a strongly-coupled medium strongly coupled to the high-p T probe. 
There are many reasons to believe that strong-coupling dynamics dominant the physics of the 
medium, and in particular, that AdS/CFT techniques provide valuable insight into these pro- 
cesses [3, 8]. For instance, running coupling calculations suggest that at T ~ 250 MeV — a not 



^t's worth noting for the experts that this ordering is not strict in pQCD: because the formation time has mass 
dependence [4] there can be violations in the ordering in certain regions of phase space. 
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Figure 1: (a) The power hq(pt) + 1 that best approximates the production spectra for various partonic species Q [31] and 
(b) Raa(pt) for pions, charm mesons, and bottom mesons at 0-20% centrality collisions [37]. Both plots use y/s^N = 2.76 
ATeV. 



unreasonable placeholder for the QGP temperature — g ~ 2 and X - g 2 N c ~ 12 » 1; it's worth 
noting that in phenomenological applications r is never large compared to Aqcd- Also, for 
T > r c , lattice calculations nontrivially deviate from the Stefan-Boltzmann limit, and in such a 
way that is reasonably well described using AdS/CFT [8]. Finally, the viscosity to entropy ratio 
extracted from hydrodynamics calculations [9] suggest rj/s ~ 1/4-n, which is readily explained 
by AdS/CFT [10]. 

There have been calculations of the energy loss of both light and heavy quarks using the 
AdS/CFT correspondence. For heavy flavor, the energy loss is a drag, dp T /dt = -fipr, where 
p = 7i A 1/2 T 2 /2Mq [11, 12]; this is similar to weakly-coupled energy loss in the Bethe-Heitler 
regime, but very different from the predictions of pQCD in the deep LPM region where dpr /dt ~ 
-LT 3 \x\(p t IMq) [13]. Comparison between AdS/CFT calculations and data are difficult be- 
cause there is no unique mapping from the parameters of QCD to those of N = 4 SYM and 
AdSs x^ 5 . Nevertheless, comparing over reasonable assumptions for parameter values in AdS/ 
CFT yields a quantitative agreement between theoretical predictions and data for non-photonic 
electron suppression at RHIC [5, 14, 15]. But we see again the power of comparing theoretical 
calculations to a wide range of data when we attempt to simultaneously describe the suppres- 
sion of heavy flavor at LHC. Keeping all parameters fixed and only changing the temperature 
of the medium (which we do according to the measured y/s dependence of the multiplicity) we 
can calculate zero parameter predictions for LHC, shown in Fig. 2. While the B meson suppres- 
sion is currently consistent with data within the large experimental and theoretical uncertainties, 
the AdS/CFT calculations significantly overpredict the suppression of D mesons. Before strong 
conclusions are drawn, however, it turns out that momentum fluctuations [16, 17] (especially 
longitudinal) — whose importance should only affect momenta parametrically large compared to 
the momenta at which the formalism breaks down, and were neglected in these calculations — 
likely play a significant role numerically. 

One wants not only to simultaneously compare AdS/CFT predictions to data as a function of 
aJs but also as a function of parton species. Unfortunately the theory of light flavor energy loss 
[21-23] is less well understood in AdS/CFT than for heavy quarks. Added difficulties arise in the 
light sector due to the lack of 1) an analytic solution for falling string configurations and 2) a good 
working definition for the energy lost by the probe (in principle one can exactly compute T^ v for 
the plasma and thus the energy lost by the probe, but this is an extremely difficult problem both 
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Figure 2: Comparison of (a) D [18] and (b) B [19] meson R A a(pt) for 0-20% centrality collisions at LHC to AdS/CFT 
heavy quark drag predictions constrained by RHIC data [5, 15, 20]. 

in terms of the analytics and the numerics). Preliminary estimates, however, suggest that light 
flavor energy loss is also overpredicted by AdS/CFT: the thermalization time for light quarks in 
the medium is of the order of 3 fm; even when the ID Hubble flow of the QGP is included the 
thermalization time is only increased to about 4 fm [24]. 

In addition to checking the flavor dependence, these AdS/CFT energy loss calculations may 
also be tested by changing collision species and looking at the suppression of heavy flavor at 
forward rapidities [25] (at mid-rapidity, open heavy flavor production can also test saturation 
effects [26]). Other information on the use of AdS/CFT in heavy ion collisions discussed during 
the conference can be found at [23, 27, 28]. 



2.2. Weakly-coupled Medium Weakly Coupled to a Probe 

The assumption of the dominance of weakly-coupled dynamics in heavy ion collisions is also 
not unreasonable. For T ~ 250 MeV, a s (2nT) = 0.3. Also, multi-loop thermal field theory is in 
good agreement with lattice data for thermodynamic properties of QCD at a few times T c , albeit 
with large uncertainties [29]. Finally, numerically daunting parton cascade calculations show 
that including 2 — > 3 channels yields rj/s ~ few/4^ [30]. 

Continuing with the assumption of a weakly-coupled plasma coupling weakly to a probe, the 
medium is described by two scales: the Debye screening length, given in terms of the Debye 
mass fi ~ gT, and a mean free path for gluons, A 8 mfp ~ l/g 2 T (see [31] and references therein). 
When evaluated at temperature scales relevant for RHIC and LHC and with all the numerical 
coefficients, one finds 2 an ordering of scales 1 Iji «c A g mfp <?c L, where again L is the pathlength 
travelled by the parent parton and is on the scale of the radius of the nucleus, L ~ R A . Since 
l/yu <?c ^ g m f p , high-/? r particles scatter off of well defined, separated medium quasi-particles. It 
is important to note however that for heavy quarks L/A m f p ~ 4 (and similarly even for gluons), 
and therefore energy loss models that assume a large number of collisions (and that thus the 
central limit theorem holds), such as those using Langevin or rates methods, likely require large 
corrections. 

In pQCD with its quasi-particle picture one can distinguish between two types of energy 
loss: elastic and inelastic, otherwise known as collisional and radiative, respectively. There is 
a long history of pQCD-based elastic energy loss calculations (see [5] and references therein); 



2 For Trhic = 350 MeV (T LHC = 450 MeV), 1/// * 0.3 (0.2) fm and A 8 mfp * 0.8 (0.7) fm. 
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see [2] for a review of pQCD-based radiative energy loss calculations. Leading order estimates 
of the size of elastic energy loss yield dpf/dt ~ -T 2 \n(p T /Mq). Naively, at asymptotically 
large energies intuition based on classical electromagnetism leads to the conclusion that AE e i «: 
AE m j, but this is based on a Bethe-Heitler estimate of radiative energy loss in which subsequent 
collisions with medium particles yield incoherently summable emissions. However, there is one 
more important scale to qualitatively understand radiative energy loss, the formation time, Tf orm , 
which characterizes the distance required for an emitted gluon to be resolved independently 
from the emitting parton. There is a large uncertainty in the size of Tf orm but for emissions of 
large energy gluons in QGP, Tf orm » ^mfp 9 anc * a sm §^ e g mon emission is thus produced from 
coherent scatterings off of multiple in-medium quasi-particles. This reduction in the amount 
of emitted radiation is known as the Landau-Pomeranchuk-Migdal, or LPM, effect [2]. In the 
LPM limit AE ra d oc -L T 3 Iti(pt/Mq). With this reduction in radiative energy loss it is possible 
for elastic energy loss to be important even at asymptotically high energies. We will be using 
the WHDG model of convolved radiative and elastic energy loss [31] for explicit comparison to 
experimental data; in this calculation which uses thermal field theory methods for computing the 
elastic energy loss and the DGLV derivation for the radiative, the elastic energy loss remains a 
significant contributor to total energy loss even for 250 GeV/c partons at LHC; see Fig. 3. 




E (GeV) E (GeV) 
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Figure 3: Elastic and radiative energy loss for gluons and light, charm, and bottom quarks travelling a distance L = 5 fm 
through QGP at (a) RHIC and (b) LHC temperatures [5]. 

We would like to compare the WHDG model to as many observables as possible. Using the 
thermal field theoretic methods to relate p and A 8 mfp to temperature, assuming that the temperature 
profile is proportional to the Glauber participant density, and that all couplings are approximately 
fixed at a s = 0.3 there is only one free parameter in the theory, the proportionality constant relat- 
ing the observed multiplicity to the entropy of the plasma. The PHENIX experiment rigorously 
extracted the best fit value of this parameter and its uncertainty, the rapidity density of gluons 
dN g /dy = 1400^ 5 [32], by comparing to their R n AA (p T ) measurement in most central = 200 
AGeV collisions. Before immediately comparing to the multitude of data from RHIC and LHC, 
it is worth noting that the lack of precision and accuracy inherent in pQCD due to both the com- 
plicated nature of the theory and the relatively large size of its coupling constant: even NLO 
calculations of production rates in hadronic collisions tend to be correct only within a factor of 
2 of the data [33, 34]. With this in mind, the agreement between the LO WHDG energy loss 
theory and data shown in Fig. 4 over a range of centralities, collision energies, measurements, 
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and flavors is surprisingly good. 
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Figure 4: Constrained zero parameter WHDG predictions compared to data for (a) V2(N part ) at RHIC [31, 35], (b) 0-5% 
centrality Raa(pt) for light flavors at LHC [34, 36], (c) R% A (pr) at 0-20% centrality at LHC [18, 37], (d) R n ^(N pa rt) at 
RHIC [31, 35], (e) v 2 (p T ) at LHC for light flavors at 40-50% centrality [37, 38], and (f) R aa (pt) at 0-20% centrality at 
LHC [19, 20]. 

There are a number of directions in which these perturbative calculations can be improved, 
many of which were discussed at the conference. For instance one might try to model the energy 
loss using a parton cascade, which trades a better treatment of multiple gluon emission for a less 
accurate treatment of the quantum mechanical formation time effects [39]. Or one might attempt 
a NLO ansatz for running coupling along with a more careful treatment of production spectra 
and time evolution [40]. For heavy flavors, some authors have attempted to include additional 
energy loss channels such as in-medium fragmentation [41] or non-perturbative cross sections 
[42]. Other heavy flavor calculations shown at this conference include [43, 44]. However, one 
should take care with calculations that include uncontrolled physics, especially ones that are 
in principle uncontrollable; it is not clear what information can be learned when unconstrained 
processes are involved. These include calculations in which parameters are dialled in to values 
that are not limited by data or are inconsistent with the assumptions underlying the model or when 
important energy loss channels are neglected (i.e. an energy loss model includes only radiative 
or elastic loss). It's also worth emphasizing that assuming a perturbative picture necessarily 
implies that the in-medium induced radiative energy loss is a result of the interference between 
the collinearly and infrared divergent radiation associated with the production of a bare color 
charge and the radiation that is created by stimulated emission from collisions of this bare color 
charge with colored objects in the QGP; calculations that do not incorporate this interference 
physics necessarily have the wrong energy and length dependencies. 



2.3. Direct Comparison of the Pictures 

Although the high-p T physics evidence for a weakly-coupled plasma weakly coupled to a 
probe is strong and there are significant signs of disagreement between the predictions of a 
strongly-coupled plasma strongly coupled to a probe it is worth considering a measurement that 
shows a qualitative difference between the two pictures. One may emphasize the different mass 
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Figure 5: Comparison of the double ratio of D meson to B meson Raa(pt) for 0-20% centrality collisions at LHC using 
the pQCD-based WHDG energy loss model [31, 36] and a model based on AdS/CFT drag energy loss [15, 20]. 

and momentum dependencies of the pQCD and AdS/CFT results by considering the double ratio 
of D to B meson Raa as seen in Fig. 5. While the leading order AdS/CFT results are applicable 
(up to a speed limits indicated on the graph) the mass dependence of the energy loss remains; on 
the other hand the mass dependence drops out for the perturbative results at asymptotically large 
momenta. 



3. Summary 

We seek a coherent, consistent picture of the physics of QGP in our quest to understand its 
properties. The comparison of energy loss calculations to data provide a direct probe of the 
relevant degrees of freedom in a QGP and how this physics interacts with high-p T particles.The 
flavor dependence is especially useful for discriminating between QGP assumptions as their 
predictions vary significantly as a function of high-p T parton masses. As demonstrated by the 
participants at this conference, heavy flavor phenomenology is rich, varied, and vigorous. 

As we have known since the first anisotropy and heavy flavor results from RHIC, a simultane- 
ous description of multiple observables related to energy loss physics is very hard to achieve. In 
particular, despite successes at RHIC, predictions for LHC based on the strong coupling physics 
of AdS/CFT do not appear to describe the data, although possibly important physics was ne- 
glected. However, LO pQCD results give a rather good qualitative description of a suite of 
observables including pion and heavy flavor suppression and anisotropy from RHIC to LHC. 
Should we find that these tentative conclusions hold, it becomes a very interesting question of 
how the strong-coupling low-p T physics of the QGP medium as implied by hydrodynamics com- 
parisons to data turn over to weak-coupling physics of the QGP medium as implied by the energy 
loss calculations. 
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